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I , I ■ The generalized Lin- Taylor model defined on the hexagonal lattice is used to investigate the phase separation in an asymmetric binary 
liquid mixture consisting of large A (hexagons) and small B (triangles) particles. By considering interaction energies between A-A 
and A-B pairs of particles that occupy nearest-neighbour cells of the hexagonal lattice, we have derived an exact solution for the 
considered model system having established a mapping correspondence with the two-dimensional Ising model on its dual triangular 
lattice. Altogether, six different types of coexistence curves including those with reentrant miscibility regions (i.e. closed-loop coexistence 
curves) were found in dependence on the relative strength between both coupling constants. 
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^ ■ 

' 1 Introduction 

', Phase equilibria of binary liquid mixtures are subject of immense research interest since closed loops 
of immiscibility were experimentally observed in the nicotine- water mixture (Q) and later onMn several 
^ aqueous solutions of the T-, 2-, 3- and 4-methylpiperidine yj l^)' /9-picoline 

i; 0; 0), 2,4- 

^ , and 2,6-dimethylpyridine 0), the glycerol-guaiacol mixture and also in many polymeric blends 
Cn| ' (|loh- It is noteworthy that an existence of the closed- loop coexistence curves with both upper (Ty) as 
: well as lower (Tl) critical solution temperatures evidently contradicts with intuitive expectations based on 

thermodynamical description of demixing, which predicts, on the contrary, a bell-shaped coexistence curve 

terminating at a single T\j instead of having an island of immiscibility. Early explanation of this remarkable 
and rather unexpected phenomenon has been suggested by Hirschfelder, Stevenson and Eyring (|llh who 
associated an appearance of the reentrant miscibility with a presence of highly orientation-dependent 
^ forces, such as hydrogen bonding, which becomes rather inefficient above a certain critical temperature. In 
_ addition, if directional bonding occurs between like as well as unlike particles and the sum of interaction 
energies for pairs of like particles is simultaneously greater than for pairs of unlike particles, the binary 
, mixture may even phase separate at low temperatures to yield an usual bell-shaped coexistence curve 
with the low-temperature Ty below an additional closed loop. Although one meets with such an intriguing 
E> ' situation rather rarely because freezing transitions usually camouflage the lowest Tu, the aqueous solution 
^ ' of 2-butanol provides a strikiiig example of the binary mixture with the closed loop above the bell-shaped 
' coexistence curve 

. 5^ 1 So far, different lattice models of liquid mixtures jiave been proposed and dealt with in order to bring 
an insight into a reentrant miscibility phenomenon (10). Even though all lattice models provide somewhat 
oversimplified picture of the liquid mixtures, some of the lattice models have already proved their worth due 
to their capability to explain the reentrant miscibility of hydrogen-bonded liquid mixtures ^16i : .17. : .lK ). Most 
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of theoretical predictions for a closed-loop formation were indeed based on the lattice-statistical models 
naturally describing the directional character of highly orientationally dependent forces and moreover, 
the lattice models can be rather easily treated using relatively precise Bethe-Guggenheim quasichemical 
approximation or Migdal-Kadanoff approach formulated on the basis of renormalization group 

theory (|22.: .2,'l: .24 : .25. : .26,). Besides, there even exist few exactly solvable lattice models, such as Frenkel- 
Louis model (1321 ) or Lin- Taylor model (LTM) which give qualitatively correct results for a phase 

separation without being affected by any approximation. It is noteworthy, however, that more complete 
understanding of the reentrant miscibility has been achieved by exploring off-lattice continuum models, 
which can be treated within the perturbation theory of Wertheim (22.; 28; .2^ or the Gibbs-ensemble 
Monte-Carlo simulation technique (jsol: Isih. 

In the present article, we shall treat one particular example of LTM recently remarkably generalized 
by Romero-Enrique and co-workers (jssl: Is^: IsTi: .38.: .39.1 . Among other matters, the generalized LTM has 
proved its usefulness in elucidating the reentrant miscibility as it has provided a plausible explanation of 
its microscopic origin. Even if the procedure worked out previously by Romero-Enrique et al. (|36l ) is rather 
general and contains our version of the model only as a special case, the numerical results reported on 
previously were mostly restricted to LTM on a square lattice, which is a rather exceptional case of the 
self-dual lattice. In this respect, we shall therefore examine here the phase separation in the binary mixture 
of large hexagonal particles and small triangular ones, which are allowed to occupy basic unit cells of the 
hexagonal lattice. As we shall see hereafter, this model also exhibits the same substantial diversity of phase 
diagrams as remarked by Romero-Enrique and co-workers (36) and all available scenarios for coexistence 
curves including those with a closed loop of immiscibility will be indeed confirmed. 

The organization of this paper is as follows. The foundations of the generalized LTM on the hexagonal 
lattice and basic steps of the exact method are given in the section |21 Section deals with the most inter- 
esting numerical results obtained for the phase diagrams and coexistence curves. Finally, some conclusions 
are mentioned in the section |3 



2 Generalized Lin- Taylor model (LTM) 

Let us begin by recalling foundations of the generalized LTM on the hexagonal lattice. Assuming that the 
binary mixture consists of two kind of species, the large hexagonal particles A and the small triangular 
ones B, LTM on the hexagonal lattice can be defined through the following rules: (i) each hexagonal unit 
cell consists of six triangular sub-units; (ii) none of A-A and A-B pairs can occupy the same hexagon; (iii) 
the small particles B can occupy the cells left empty by the large particles A; (iv) each triangular sub- 
unit can be occupied by at most one particle B; (v) pairwise interactions saa and sab exist between the 
nearest- neighbouring A-A and A-B pairs sharing a common edge (Fig. 1). Next, we shall refer to /xa and 
jiB as to the chemical potentials relative to the particles A and B, respectively. For further convenience, 
let us also define fugacities of both the particles by 

ZA,B = exp(/3/iA,B), (1) 

where (3 = \/{k-QT), k-Q is Boltzmann's constant and T absolute temperature. The grand-canonical parti- 
tion function of this binary liquid mixture can be written in the form of a sum taken over all occupation 
configurations {nj} available for a given system 

H = ^ exp(^/?/XA ^^ni- (3eAA ^ '^B{{ni}). (2) 

Above, rij = 1,0 labels occupation number of the basic unit cell (hexagon), whereas nj = 1 if ith hexagon 
is occupied by the large particle A, otherwise iii = 0. The summation Yli runs over all basic unit cells of 
the hexagonal lattice and the other one j) is performed over all pairs of nearest-neighbour unit cells. If 
the configuration of large hexagonal particles A is fixed, the grand-canonical partition function SB({^^^}) 
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of the small particles B can be expressed as follows 

HB({na) = (1 + zb)*'^(1 + ^Be-^^*-)^^ (3) 

Here, Mp (Mi) denotes the total number of triangular sub-units available for occupancy to the small 
particles B, which have (not) a common edge with the large particles A. When taking into account the 
lattice geometry, both the numbers can easily be evaluated in terms of the occupation numbers 

Mf = 6^ni-2^ninj, (4) 
Mi = 6N - 12^^711 + (5) 

where labels the total number of the lattice cells (hexagons). At this stage, the grand-canonical partition 
function of the binary mixture can be mapped onto an one-component lattice-gas model for the species A 
by substituting (jHJ, (HJ) and © into the relation Q and thus, performing a summation over the occupation 
variables relative to the species B 

H = (1 + ZBf^ Yl exp{[-/3eAA - 21n(l + ZBe"^"*") + 21n(l + zb)] X^^^^j} 

X exp{[/3/iA + 61n(l + ZBe-^^*^)-121n(l + ZB)] J^ni}. (6) 

i 

The equality (jHl) proves that the grand-canonical partition function can be calculated from the equiva- 
lent lattice gas model with the renormalized pair interaction and shifted chemical potential. Next, the 
transformation between the occupation number (n^ = 1,0) and Ising spin variable (sj = ±1) 

n, = (1 + s,)/2, (7) 

establishes a mapping relationship between the model under investigation and the spin-^ Ising model 
on the dual triangular lattice. As a consequence of this, the grand-canonical partition function can be 
expressed through the partition function of the Ising model on the triangular lattice 

E = [4e-3/^^**(l + zb)^(1 + ZBe-^^*-)6]f Zt,iang(i2,i^). (8) 

Note that the mapping parameters R and H stand for the effective coupling constant and effective magnetic 
field in the associated spin-^ Ising model on the triangular lattice 



fi=-lne-^"^ \ ^ ,J (9) 

1 /_~3/3e>4 



2 V (1 + zb)V ^ ^ 

It is worthwhile to remark that only the requirement of zero effective field {H = 0) in the equivalent 
Ising model ensures a possible existence of the phase separation in the considered LTM. As a result, the 
necessary (but not sufficient) condition allocating the coexistence region between two phases of diff'erent 
composition relates two fugacities of the large and small particles through za = e^'^'''^*(l -|- zb)^- 

Notice that the mapping relation ^ represents a central result of our calculation as it formally completes 
an exact solution for the grand-canonical partition function H with regard to the known exact result of the 
partition function of the triangular Ising model (j4ol : 41). Besides, this relationship can readily be utilized 
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for calculating the number density of the large {n^) and small (ne) particles by relating them to the 
magnetization per site ( mtriang-) and the internal energy per site (ntriang) of the corresponding Ising model 
on the triangular lattice (|4fl : 4ll : 43 ). After straightforward calculation one actually finds 



TiA = ^ ^ ^*"^"S ^ rag = ^(C^ + C2) - SCimtriang " ^(Cl - C2)ntriang, (H) 

where the coefficients Ci and C2 are defined as follows 

Ci = -^, C2 = -^^5 . (12) 

The composition of the binary mixture is subsequently unambiguously determined by the molar fraction 
of the large particles A defined by means of 

Xa = (13) 
Finally, the pressure of the binary mixture can be calculated from the grand potential by the use of 

where a labels the area of basic unit cell (hexagon) . For simplicity, we shall use throughout the rest of this 
paper the renormalized quantity p = p*Na as a measure of the pressure. 



2.1 Coexistence region 

In order to describe the phase separation into two phases of different composition, let us firstly express the 
fugacity and pressure as a function of temperature, the interaction parameters eaa, £ab and the effective 
coupling R. According to the equation ® , the fugacity of small particles B can be calculated from 

g2iJ _ g-/3£AA/2 

^ e-/3£AA/2 _ e2i?-/3£AB ■ ^^^"^ 
Furthermore, substituting (|15j) and ^ into the equation H14|) gives the following expression for pressure 

\ _ g-/3£AE 



P 



r r 1 — e~'^^^*^ 1 1 1 

kBT[-3R + 6\n[ ^.2i;-/3.AA/2 _ ^-Pe.^ J + jV ^" ^^"""^ I " ^^^^ 



Experimental conditions under which coexistence curves are usually observed meet the requirement of 
constant pressure. In this respect, the equation (jlOj) can serve for determining the effective coupling 
parameter R by selecting T, eaa and eab at fixed value of pressure. A subsequent substitution of the 
effective coupling R into the equation (|15() then enables a simple calculation of the fugacities (the fugacity 
of large particles is connected to the one of small particles due to an unconditional validity of the zero-field 
condition). With all this in mind, the composition of binary liquid mixture can be consequently calculated 
using the set of equations (fTT]) - ((T^ . 



2.2 Critical behaviour 

According to the mapping relation Q, the binary mixture becomes critical if, and only if, the effective 
coupling parameter of the associated Ising model approaches its critical value. Owing to this fact, the 
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critical points can be found from this simple equation 



0.055627 + 6 In 



r ^3(1 -e-^-^*«) 



} 



(17) 



B 



where Pc = l/(^B^c)) Pc and Tc label the critical pressure and critical temperature. The molar fraction of 
large particles at a critical point can be consecutively simplified to 



3 Results and Discussion 

In this part, we shall briefly discuss the most interesting results obtained for the phase diagrams and 
coexistence curves of the generalized LTM on the hexagonal lattice. It is worthy to mention that the 
displayed coexistence curves are obtained by cutting concentration-temperature slices at fixed pressure, i.e. 
a situation which fully corresponds experimental conditions by performing studies of the phase separation. 
It directly follows from the equation (|15|) that the phase separation into A-rich and B-rich phases might in 
principle appear either if saa > 0, eab/^aa > 0.5, or ^aa < (arbitrary £ab)- In what follows, we shall 
treat those special cases in several sub-sections. 

3.1 Repulsive interaction £aa > 

If the pair interaction between the particles A is repulsive, the coexistence then possibly occurs just as 
the interaction between the A-B pairs is likewise repulsive and simultaneously, the repulsion energy eab is 
stronger than a half of the repulsion energy eaa, i-e. eAe/^AA > 0.5. The phase diagram for this particular 
case is depicted in Fig. 2a in the form of Tc — Pc dependence. The region inherent to the coexistence 
of the A-rich and B-rich phases can be located above the displayed curves, while below them both the 
components become perfectly miscible. It is quite obvious from Fig. 2a that the phase separation occurs 
only above a certain boundary pressure Pb/^AA = 3, which is needed for overcoming the repulsive force 
between the pairs of particles A. On the other hand, the mixture becomes perfectly miscible independently 
of the ratio sab/^aa below this pressure value. It can be also readily understood from Fig. 2a that the 
stronger the repulsion energy between the A-B pairs, which means, the stronger the ratio eAe/eAAi the 
higher the critical temperature at which the phase coexistence disappears at a given pressure. Next, the 
lines depicted in Fig. 3a illustrate changes of the critical concentration along the critical lines from Fig. 
2a. As one can see, the critical concentration achieves = ^ regardless of eab/^aa when approaching 
the lowest pressure pb at which the phase separation appears. For completeness, three typical coexistence 
curves are shown in Fig. 4a for eab/^ka = 1 fixed and several values of pressure. As it can be clearly seen, 
one finds the typical bell-shaped coexistence curves irrespective of the pressure strength. 

3.2 Attractive interaction eaa < 

Now, let us examine the case when the pairwise interaction between the particles A is attractive. In such 
a case one encounters much more complex phase diagrams and also much richer critical behaviour. For 
better orientation, we have divided our discussion into several parts dealing with some typical cases. 

3.2.1 £ab > 0. In this case, the repulsive interaction between the A-B pairs favours the phase separation 
into A-rich and B-rich phases and as a consequence of this, one finds a quite similar phase diagram and 
coexistence curves as discussed in above (compare Fig. 2b with 2a). The most significant difference consists 
in an appearance of a special critical end point with coordinates [k-oT* /\eAA\^ P*c/\^ A a\\ = [0.910,0.051], 




(18) 
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which is labelled by a star symbol in Figs. 2b)-d). While for pressures greater than p* the critical tem- 
perature monotonically increases as pressure rises, below the pressure there is a coexistence but no 
criticality. This observation would suggest that the critical end point terminates at a liquid-vapour coexis- 
tence line of the pure A component. Actually, the mixture becomes pure A before it turns critical for any 
pressure lower than p*. This fact can be clearly seen also from Fig. 3b, where the critical concentration 
is plotted against the critical temperature and all lines start from the same critical temperature of the 
pure A component (X^ = 1). For illustration, we depict in Fig. 4b three typical coexistence curves that 
obviously exhibit much more pronounced asymmetry than those discussed formerly for the case with the 
repulsive interaction ^aa- Apparently, this asymmetry is the more evident, the lower and the closer is 
pressure selected to its critical value p* , since the initially B-rich phase looses at low pressures much more 
rapidly its B component than the A-rich phase is enriched by the particles B. 



3.2.2 —0.5 < £ab/|£aa| < 0. Contrary to the aforedescribed behaviour, it is easy to understand from 
Fig. 2c that a weak attractive force between the A-B pairs leads to a monotonous decrease of the critical 
temperature when increasing the pressure strength. The standard bell-shaped coexistence curves, which 
are plotted in Fig. 4c for one particular ratio eAs/kAAl = —0.3 and several values of pressure, provide 
a strong support for this statement. According to these plots and also dependences drawn in Fig. 3c, a 
suppression of the critical concentration in response to a pressure strengthening is observed due to a 
change of the character of the interaction energy ^ab- This rather easily understandable behaviour survives 
unless the ratio between both coupling constants does not reach the value eAB/kAA| = —0.431. Within 
the interval —0.5 < eAe/kAAl < —0.431, however, an outstanding non-monotonous dependence of the 
critical temperature on pressure can be detected. Assuming for instance that eAB/kAA| = —0.45 is kept 
constant, the —pc dependence can be characterized by one local minimum [kBT"^"^/\eAA\,pT"^/\^AA\] = 
[0.146,1.436] and one local maximum [kBT^''''/\eAA\,pT''''/\£AA\] = [0.353,1.643]. While below pf"" or 
above p™^^ the usual bell-shaped coexistence curves with single T\j should be expected, the coexistence 
curves with three consecutive critical points (Tu, Tl, and T\j) and reentrant miscibility should emerge 
for any pressure from inside the interval bounded by p™™ and p^^^. Fig. 4d illustrates such example of a 
closed loop above a bell-shaped coexistence curve with in total three critical points obtained at pressure 
p/\£aa\ = 1-5, as well as, the usual bell-shaped coexistence curve to emerge whenp/|eAA| = 0.5. If pressure 
is selected sufficiently close but still below p™™, the 'hour-glass' coexistence curve can be even detected; 
temperature induces after initial increase of mutual solubility its decrement until a repeated increase of 
solubility near Tu re-appears, as it is clearly depicted for the particular case p/\eAA\ = 1-0. It is worthwhile 
to remark that the aforementioned behaviour is completely consistent with experimental observations of 
the pressure effect on the coexistence curves of the aqueous solution of 2-butanol (0; 3 ish - 

Even more involved situation emerges if pressure is selected directly equal to p™™ or p™*^^. When the 
upper bound p'^^^ is selected, then Tu and Tl incident to an island of immiscibility coalesce at so-called 
double critical points which are shown as circles. On the other hand, the low-temperature Tu of bell- 
shaped curve merges together with Tl of the closed loop to yield a critical double point (which is marked 
as diamond) by selecting the lower bound p™™ for pressure. Notice that the double critical points as well as 
the critical double point can be characterized by a doubling of the critical exponents as it has been proved 
previously (js^)- However, the most peculiar critical point [ksTc^ /\£aa\,p'c /\£aa\] = [0.265,1.862] appears 
by considering following value for the ratio between both coupling constants eAB/kAA| = —0.431. In this 
case all three critical temperatures coalesce simultaneously at so-called critical inflection point, which can 
be characterized by a tripling of the critical exponents (|37l). 



3.2.3 eaa < 0, £ab/|£aa| < —0.5. Finally, we shall discuss the coexistence phenomenon for the case 
when the attractive force between the A-B pairs is stronger than a half of the attractive force between 
the A-A pairs, which means, eAB/|eAA| < —0.5. Under these circumstances, the coexistence region can 
be located below the curves displayed in Fig. 2d. It is quite obvious from this figure that there appears a 
closed loop of immiscibility whenever pressure is selected above p* but below the value p^^^ corresponding 
to the double critical point. Namely, the miscibility gap A = Tu — Tl gradually decreases upon pressure 
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strengthening until it vanishes at the double critical points where Tu and Tl merge together. This scenario 
can be clearly seen from Fig. 4e, where the closed-loop coexistence curves are plotted for three different 
pressures and the ratio eAB/|eAA| = —1- It is quite apparent from this figure that the stronger the 
pressure, the smaller the miscibility gap (coexistence region). Note furthermore that there even exists 
another possible scenario to have a coexistence without criticality at Tu when selecting pressure below 
p*/|eAA| = 0.051. In such a case, the temperature induces a phase separation at Tl, nevertheless, both 
separated phases become pure A before the mixture turns critical and hence, one observes two pure phases 
of particles A in coexistence that merely differ in their densities. This observation can be interpreted as 
a liquid-vapour phase separation of the pure A component, which takes place because the component B 
vapourises prior to achieving the critical temperature when pressure is selected below its critical value p*. 

4 Conclusion 

The present article deals with the investigation of phase separation in an asymmetric binary liquid mixture 
of large hexagonal and small triangular particles described by means of the generalized LTM on the 
hexagonal lattice. Despite its simplicity and a certain oversimplification, this model is relevant as it provides 
deeper understanding of the phase separation phenomenon without making any further approximation to 
the results obtained. In addition, this model proves an existence of the closed loops of immiscibility 
occurring under certain conditions either separately, or above the standard bell-shaped coexistence curves. 
The closed-loop coexistence curve bounded by two critical points Tu and Tl indicates such a kind of 
the reentrant miscibility, where Tl determines an upper critical temperature under which two components 
become repeatedly perfectly miscible. On the other hand, the reentrant miscibility can also be found in the 
closed loop plus bell-shaped coexistence curve with in total three critical points; the region of reentrant 
miscibility then occurs in between Tl of the closed loop and Tu of the bell-shaped coexistence curve. 
Altogether, six possible scenarios for the phase separation were illustrated with the help of exact results 
for the phase diagrams and coexistence curves: the standard bell-shaped curve with Tu, the bell-shaped 
curve without Tu, the bell-shaped curve plus closed loop with two Tu and one Tl, the closed- loop curve 
with Tu and Tl, the 'hour-glass' curve with Tu, and the closed loop with Tl but without Tu- 

The main objective of the present work was to provide an eventual confirmation of the aforementioned 
coexistence scenarios, which were originally envisaged by Romero-Enrique and co-workers after introducing 
and exploring the generalized LTM (|36h. It should be pointed out, however, that the majority of numerical 
results reported on hitherto were mostly restricted to the particular case of the generalized LTM on the 
square lattice. From this point of view, the investigation of LTM on the hexagonal lattice is of a particular 
importance, because this model even has more obvious relevance to the phase separation of real binary 
liquid mixtures which might consist of molecules with a hexagonal symmetry (like benzene, cyclohexane 
and a large number of their structural and heterocyclic derivatives) and smaller non-linear molecules of 
the solvent (like acetone, isopropanol, ethyleneoxide, ethers, dimethylsulf oxide, etc.). Our next effort is to 
provide a further extension to the model under consideration to account also for multiparticle interactions 
and to elucidate a role of the size of solvent on the mutual solubility. 
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Figure Captions: 

Fig. 1 One among possible configurations of particles within LTM on the hexagonal lattice. Large grey hexagons 
denote lattice positions of the particles A, small black triangles positions of the particles B. The inter- 
actions are considered only between A-A and A-B pairs of nearest neighbours, which share a common 
edge (point contacts are neglected). 

Fig. 2 Phase diagrams in the form of Tc — Pc dependence for four different cases considered: a) ^aa > 0, 
^AB/eAA > 0.5; b) £aa < 0, eAB/|eAA| > 0; c) £aa < 0, -0.5 < sab/IsaaI < 0; d) eaa < 0, 
£ab/|£aa| < —0.5. The star symbol denotes a special critical end point terminating on a liquid-vapour 
coexistence line of the pure A component, the circles and diamond label the double critical points and 
critical double point, respectively. For clarity, the points at which coexistence of two different pure A 



February 2, 2008 



7:52 Molecular Physics strecka 



REFERENCES 9 

phases disappears prior to achieving the critical point is shown as a broken hne only in Fig. 2d) for the 
particular ratio eAB/|eAA| = — 1- 

Fig. 3 The changes of critical concentration Xj^ along the critical lines displayed in Fig. 1. The star symbol 
denotes the critical end point terminating on a liquid- vapour coexistence line of the pure A component. 

Fig. 4 Some typical examples of the coexistence curves displayed in the form of composition vs. temperature 
dependence for: a) £aa > 0, eAA/^AB = 1, b) eaa < 0, £:aa/|£ab| = 1, c) Eaa < 0, eAA/|£AB| = -0.3, 
d) £aa < 0, eAA/|eAB| = —0.45, e)-f) eaa < 0, eAA/kAB| = ^1, and several values of pressure. Open 
circles determine a position of the critical points, broken lines display their pressure changes. 



